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SUMMARY 


// 

Rate  gyroscope  biases  play  an  important  role  in  flight 
tests  requiring  flight  path  reconstruction,  a  method  often 
used  in  aircraft  parameter  estimation.  The  biases  can  drift 
with  time,  be  affected  by  system  power  up  and,  ideally, 
should  be  calibrated  before  e^ch  flight  to  maintain  optimum 
performance.  This  report  details  a  numerical  method  to 
determine  the  biases  of  high  quality  flight  test  rate  gyro¬ 
scopes  immediately  prior  to  a  flight.  The  accuracy  of  the 
method  is  such  that  the  earth  rate  is  clearly  sensed  and 
accounted  for,  a  variable  rarely  considered  in  flight  test¬ 
ing.  The  method  requires  minimal  time  and  no  calibration 
hardware .  ,■/ 


RESUME 

Les  erreurs  systematiques  des  gyrometres  ont  une  grande 
influence  dans  les  essais  en  vol  qui  necessitent  la 
reconstitution  d'une  trajectoire  de  vol,  methode  souvent 
utilisee  pour  1’ estimation  des  parametres  d'un  aeronef.  Les 
erreurs  systematiques  peuvent  varier  avec  le  temps  et  peuvent 
etre  influencees  par  la  mise  en  marche  de  1 '  alimentation  du 
systeme;  c'est  pourquoj ,  idealement,  les  gyrometres  devraient 
etre  etalonnes  avant  chaque  vol  afin  qu'ils  puissent  offrir 
des  performances  optimales.  Le  present  rapport  decrit  une 
methode  numerique  pour  determiner  les  erreurs  systematiques 
des  gyrometres  de  precision  immediatement  avant  le  vol.  La 
precision  de  cette  methode  est  telle,  que  le  taux  de 
precession  apparente  de  la  terre  est  clairement  detecte,  et 
que  cette  variable  est  prise  en  compte  dan3  les  calculs,  ce 
qui  est  rarement  le  cas  pour  les  essais  en  vol.  La  methode 
n'exige  que  peu  de  temps,  et  ne  fait  appel  a  aucun  instrument 
d' etalonnage. 
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NOMENCLATURE 


Symbol  Definition 

a  initial  x-accelerometer  measurement  (average) 

o 

a  initial  y-accelerometer  measurement  (average) 

yo 

a  initial  z-accelerometer  measurement  (average) 

*0 

bp  roll  rate  gyro  bias 

bg  pitch  rate  gyro  bias 

br  yaw  rate  gyro  bias 

L  geocentric  latitude,  positive  in  northern  hemisphere 

p  true  angular  velocity  component  in  roll 

(positive  direction  is  right-wing  down) 

q  true  angular  velocity  component  in  pitch 

(positive  direction  is  nose  up) 

r  true  angular  velocity  component  in  yaw 

(positive  direction  is  nose  to  the  right) 

Pm  measured  angular  velocity  component  in  roll 

qm  measured  angular  velocity  component  in  pitch 

rm  measured  angular  velocity  component  in  yaw 

8  pitch  angle  (positive  nose  up) 

0  roll  angle  (positive  right  wing  down) 

♦  yaw  angle  (see  figure  2  for  reference  system) 

Qe  Earth's  rotation  rate  relative  to  inertial  space 

geocentric  longitude 
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AN  ACCURATE  NUMERICAL  TECHNIQUE  FOR  DETERMINING  FLIGHT  TEST 
RATE  GYROSCOPE  DIASES  PRIOR  TO  TAKEOFF 


1.0  INTRODUCTION 

In  the  context  of  flight  testing,  flight  path  recon¬ 
struction  refers  to  integration  of  the  inertial  measurements 
over  the  brief  time  period  of  the  manoeuvre,  generating  time 
histories  of  aircraft  attitude  and  velocity  components.  The 
typical  length  of  a  manoeuvre  is  approximately  1  minute  or 
less,  hence  the  requirements  for  measurement  accuracy  are 
substantially  less  than  that  for  inertial  navigation  where 
integration  times  of  more  than  1  hour  are  often  required 
(with  measurements  being  augmented  by  sophisticated  Kalman 
filter  integrated  navigation  algorithms  to  eliminate  dominant 
errors).  The  flight  path  reconstruction  results  are,  none¬ 
theless,  significantly  affected  by  the  instrument  biases,  and 
considerable  effort  is  often  spent  in  determining  their 
value. 

The  University  cf  Toronto  Institute  for  Aerospace  Stud¬ 
ies  flight  test  data  acquisition  system  was  obtained  and  used 
in  this  investigation.  The  temperature  controlled  inertial 
module  contains  3  Honeywell  GG1111AN02  rate  gyroscopes  and  3 
Sundstrand  QA-2000  Q-Flex  accelerometers.  This  is  considered 
very  high  quality  fliqht  test  instrumentation. 

To  maintain  peak  performance,  these  instruments  should 
be  routinely  calibrated  to  identify  drift  with  time. 
Although  impractical,  they  should  be  calibrated  after  each 
turn-on  (before  each  flight)  to  establish  the  turn-on  to 
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turn-on  bias  shift.  Furthermore,  many  data  acquisition  sys¬ 
tems  contain  analogue  to  digital  converters  that  are  subject 
to  bias  drift  with  time,  which  may  be  interpreted  as  an  in¬ 
strument  bias. 

This  report  addresses  a  practical  and  accurate  method  of 
determining  rate  gyroscope  biases  prior  to  each  flight. 


2 . 0  TEST  PROCEDURE 


Although  the  following  technique  discusses  the  setup 
within  the  laboratory,  it  cam  easily  be  extended  to  include 
implementation  within  the  aircraft  prior  to  flight.  The  main 
assumptions  are  that  the  principal  error  sources  are  biases, 
and  that  they  remain  constant  for  the  duration  of  the  experi¬ 
ment  (=  90  seconds). 

The  technique  requires  an  inertial  module  capable  of 
being  rotated  in  all  directions,  and  has  the  following  simple 
steps : 

1)  Place  the  inertial  module  on  any  3tatiorary  sur¬ 
face. 

2)  Begin  recording  data. 

3)  Manually  rotate  the  inertial  module  about  each  of 
the  three  axes  in  both  the  positive  and  negative 
directions.  Care  must  be  taken  so  as  not  to  over¬ 
range  the  instruments.  See  sample  time  histories  of 
figure  1. 

4)  Return  the  inertial  module  to  the  starting  location 
(exactly) . 

5)  Stop  recording  data. 


The  idea  is  to  select,  via  a  searching  algorithm,  the 
gyro  biases  required  to  compensate  for  the  differences  be¬ 
tween  indicated  values  of  the  starting  and  ending  angles  of 
pitch,  roll  and  yaw  are  zero.  The  search  algorithm  can,  cf 
course,  be  executed  post-flight  since  this  is  typically  when 
all  data  analysis  is  conducted.  Knowledce  of  the  initial 
pitch  angle  and  roll  angle  is  obtained  with  the  accelerometer 
readings  or  with  an  inclinometer.  The  initial  yaw  angle 
(azimuthal  direction  from  East)  is  required  to  account  for 
the  Earth's  rotation  as  will  be  demonstrated  in  section  4.1. 

3.0  NUMERICAL  PROCEDURE 

The  numerical  procedure  described  in  sections  3.1  and 
3.2  has  been  implemented  on  the  NRC  IBM  3090  mainframe  ir. 
double  precision.  Total  execution  time  for  a  90  second  lab¬ 
oratory  "manoeuvre"  (40  Hz  sampling  rate)  is  approximately  35 
CPU  seconds.  Although  the  procedure  is  iterative,  no  conver¬ 
gence  problems  were  encountered. 

3 . 1  RATE  EQUATIONS 

The  complete  equations  expressing  the  relationship  be¬ 
tween  rotation  rate  and  attitude  are  taken  from  reference  1: 


* 


rcosff  qs.4  ig 
COS0  COS0 


♦  Lcos^tan© 


♦  (Q  *  X)(sinL  ♦  cosLs intf«tan0 ) 

© 


3.1.1 


0  =  qcos0  -  rsin0  -  [Lain*  -  \ Q  +  X)cosLcos*>] 


3.1.2 


4 


=  p  +  qtan9sin0  »  rtanOccp^J 

♦  r  _Lcos*  +  (p  +  i)cogL  sir !± 
l  roafl  '  e  '  cos0 


3.1.3 


For  this  experiment,  the  rate  of  change  in  latitude  and 
rate  of  change  in  longitude  (L  and  i)  are  of  course  zero 
(flat  Earth  assumption),  hence  the  equations  can  be  simpli¬ 
fied  to: 


*  -  rco.?$.  ♦  ♦  Q  (sinL  +  cosLsin*tan0 ) 

cos0  cos9  ev 


3.1.4 


=  qcos0  -  rsin?l  ♦  Q  cosLcosij) 


3 . 1 .  S 


-  p  ♦  qtan9sin 0  ♦  rtan0cos{3  >•  QecosL  — 


3.1.6 


The  Earth's  rotation  rate  (relative  to  inc-rtial  space) 
is  15.04107  deg/hr,  reference  2.  When  the  Earth  rate  is 
accounted  for,  the  azimuthal  direction  of  the  inertial  module 
is  required.  The  reference  system  is  such  that  when  the  in¬ 
ertial  module  (or  aircraft)  is  pointing  East,  *  -  0.0  de¬ 
grees.  The  sign  convention  is  such  that  when  the  inertial 
module  is  facing  North,  ♦  =  -90.0  degrees  (see  figure  2). 

When  accounting  for  the  rate  gyroscope  biases,  the  equa¬ 
tions  become: 


<r„-b^)coa0  ain0 

i  s  £ £ -  2 — 2- -  ♦  Q  (sinL  ♦  cosLsinftanG ] 

*  co80  CO80  e' 


3.1.7 
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0  =  (q  -  b  )cos0  -  (r  -  b  )sin0  51  cosLcos*  3.1.8 

m  m  r  e 


(Pm  -  bp>  ♦  <<V  bq)tan0sin0  ♦  <rm 

♦  c  cosL  1114 


b^ ) tan0cos0 


3.1.9 


where  L  =  45.325509  degrees 


These  are  the  equations  used  in  this  study. 


3.2  INTEGRATION  AND  SEARCH  METHOD 

In  order  to  determine  the  biases,  the  above  equations 
must  be  integrated  from  the  start  of  the  laboratory  "ma¬ 
noeuvre"  to  the  end.  The  pitch,  roll  and  yaw  angles  at  the 
end  of  the  integration  are  compared  to  their  values  at  the 
beginning.  If  the  comparison  does  not  result  in  the  begin¬ 
ning  and  ending  angles  being  the  sar.  e  for  each  of  the  three 
axes,  the  three  rate  gyroscope  biases  are  appropriately  ad¬ 
justed  to  force  agreement  -  an  iterative  procedure. 

The  initial  pitch  and  roll  angles  used  during  the  inte¬ 
gration  can  be  obtained  from  the  accelerometer  readings,  as 
shown  below.  These  angles  are  not  required  with  high  accu¬ 
racy  since  they  have  minimal  influence  on  the  resulting 
biases.  For  example,  the  largest  error  in  bia3  as  a  result 
of  a  5  degree  error  in  initial  pitch  angle  was  only  0 . 5 
deg/hr  (Q-gyro),  as  determined  in  a  numerical  experiment. 
This  appealing  insensit : vity  to  initial  conditions  results 
from  the  fact  that  the  numerical  constraint  is  the  change  in 


angle  (or  lack  thereof),  not  the  absolute  value. 
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0 

o 


s  m 


-1 


0  =  sin 

o 


or : 


gcos6o 


) 


-ax  cosJ3o 

0  =  tan  (  - — - 

o  a„ 


^o  = 


tan 


-1 


3.2.1 


3.2.2 


3.2.3 


3.2.4 


The  inertial  module  remains  stationary  for  the  first  5 
seconds  of  data  recording,  to  allow  average  accelerometer 
readings  to  be  calculated  and  used  in  the  equations  above. 
The  error  introduced  by  the  accelerometer  biases  in  the  above 
equations  will  be  small,  since  integration  over  time  is  not 
involved  (i.e.  a  bias  of  1  milli-g  will  produce  an  angular 
error  of  only  3.4  arc  minutes).  An  inclinometer  can  also  be 
used  to  determine  the  inertial  modules  initial  pitch  and  roll 
angles.  The  initial  yaw  angle  must  be  known  since  the  Earth 
rate  sensed  by  the  individual  rate  gyroscopes  is  generally  a 
function  of  azimuthal  position.  Since  accounting  for  the 
Earth  rate  is  considered  a  small  correction  for  flight  test 
purposes,  a  small  error  in  the  initial  yaw  angle  is 
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unimportant.  A  numerical  experiment  has  shown  that  an  error 
in  initial  yaw  angle  (?0)  of  5  degrees  results  in  a  P-gyro 
bias  error  of  only  2  deg/hr. 

The  integration  method  used  in  this  study  was  a  5th  and 
6th  order  Runge-Kutta-Verner  pair  (DIVPRK  of  reference  3). 
The  routine  used  to  conduct  the  nonlinear  search  for  the 
biases  was  a  Levenberg-Marquardt  algorithm  (DNEQNF  of  refer¬ 
ence  3).  This  algorithm  routinely  converged,  even  with  ini¬ 
tial  bias  estimates  set  to  zero.  A  flowchart  of  the  pro¬ 
cedure  is  given  in  figure  '* . 

4.0  DISCUSSION  OF  RESULTS 

A  total  of  21  separate  "manoeuvres"  were  conducted  on 
three  separate  days  within  a  two  week  period.  The  experi¬ 
ments  were  conducted  on  separate  days  to  check  for  turn-on  to 
turn-on  repeatability  of  the  biases. 

The  first  turn-on  consisted  of  four  "manoeuvres",  the 
results  of  which  are  presented  in  figure  4.  For  each  of  the 
three  rate  gyroscopes,  the  computed  biases  are  seen  to  be 
very  repeatable. 

The  biases  of  flight  test  gyroscopes  prior  to  flight  are 
often  estimated  by  simply  recording  data  with  the  inertial 
module  stationary.  Ideally,  the  rate  gyroscopes  should  read 
zero  under  these  conditions  (excluding  the  Earth  rate  -  which 
is  normal  practice  in  flight  testing),  hence  any  non-zero 
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reading  is  considered  a  bias.  Note  that  for  the  UTIAS  12  bit 
A/D  conversion  (4096  counts)  and  ±  25  deg/sec  rnnge  lor  each 
channel,  the  finest  resolution  achievable  is  0.0122  deg/sec 
or  44  deg/hour.  Hence,  the  biases  can  only  be  determined  to 
a  resolution  of  44  deg/hour  by  using  the  above  stationary 
module  method.  The  method  presented  herein  has  essentially 
infinite  resolution  and  the  results  presented  in  figure  4  are 
within  the  44  deg/hour  resolution  of  the  stationary  module 
method  for  all  three  gyroscopes.  The  results  indicate  that 
biases  equivalent  to  2  or  3  A/D  conversion  increments 
(counts)  are  present  on  all  three  channels.  Three  candidate 
sources  of  these  relatively  large  biases  are  initial  calibra¬ 
tion  error,  rate  gyroscope  bias  shift  after  calibration  and 
A/D  converter  shift.  The  latter  is  considered  to  be  the  most 
probable. 

The  results  from  the  second  and  third  turn-ons  ar“  pre¬ 
sented  in  figures  5  and  6,  which  are  similar  to  the  results 
from  the  first  turn-on.  The  average  biases  for  the  three 
turn-ons  are: 


Q-rate  gyro 
P-rate  gyro 
R-rate  gyro 


TURN-ON  1 
(deg /hr) 

-100.9 
-  79.9 
-108.1 


TURN-ON  2 
(deg/hr) 

-  99.0 

-  75.6 
-107.6 


TURN-ON  3 
(deg/hr) 

-101.2 
-  77.4 
-109.2 


It  is  clear  that  the  variation  in  bias  from  turn-on  to 
turn-on,  for  a  relatively  short  2  week  period,  is  minimal. 


9 


The  results  for  all  21  "manoeuvres"  are  presented  in  figure 
7.  The  overall  averages  and  standard  deviations  for  the 
three  gyroscope  biases  are: 


AVERAGE 

(deg/hr) 

Q-rate  gyro  -100.2 

P-rate  gyro  -  77.1 

R-rate  gyro  -108.3 


STANDARD  DEVIATION 
(deg/hr) 

1.9 

1.9 

2.6 


A  standard  deviation  of  roughly  2  deg/hr  (0.00055  deg/ 
sec)  is  considered  excellent,  since  the  Earth  rate  (approxi¬ 
mately  10.6  deg/hr  at  45  degrees  latitude)  is  rarely  even 
considered  for  flight  test  purposes  during  standard  calibra¬ 
tions.  Results  of  this  quality  indicate  that  the  initial 
assumptions  are  valid,  that  is,  the  biases  are  the  dominant 
error  sources  and  that  they  remain  essentially  constant  for 
the  duration  of  the  manoeuvre.  The  results  also  demonstrate 
the  stability  (in  the  short  term  at  least)  of  the  UTIAS  rate 
gyroscopes  and  data  acquisition  system. 

To  achieve  results  with  2  deg/hour  standard  deviation 
requires  highly  repeatable  starting  and  ending  points.  Con¬ 
sider,  for  example,  that  if  the  true  ending  angle  in  pitch 
was  different  than  the  true  starting  angle  by  a  0.05  degree 
misplacement  of  the  inertial  module  (3  arc  minutes),  this 
alone  introduces  an  error  of  2  deg/hour  for  a  90  second 
"manoeuvre".  Hence,  care  must  be  exercised  in  the  implemen¬ 
tation  of  the  method. 
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4.1  EFFECT  OF  EARTH  RATE  ON  RESULTS 

In  order  to  demonstrate  that  the  technique  proposed 
herein  was  of  sufficient  accuracy  to  warrant  the  inclusion  of 
the  Earth  rate  in  the  calculations,  a  special  test  procedure 
was  conducted  during  the  third  turn-on. 

For  the  first  two  turn-ons,  13  "manoeuvres",  the  begin¬ 
ning  and  ending  azimuthal  positions  were  the  same,  *=-33.17 
degrees,  and  the  "manoeuvre"  time  histories  were  very  simi¬ 
lar.  Each  of  the  three  gyroscopes  sensed  a  different,  but 
repeatable,  component  of  Earth  rate  as  shown  below: 

q  =  -flecosLcos*  4.1.1 

p  =  -QecosLsin*  4.1.2 

r  =  -Q  sinL  4.1.3 

e 

These  are  the  components  of  Earth  rate  sensed  by  a  sta¬ 
tionary  inertial  module  whose  pitch  and  roll  angles  are  zero. 
Note  that  azimuthal  orientation  affects  the  Q  and  P  rate 
gyros  and  not  the  R  gyro  under  these  conditions.  For  a  par¬ 
ticular  beginning  and  ending  azimuthal  orientation  and  simi¬ 
lar  "manoeuvre"  time  histories,  the  resulting  biases  would 
have  been  repeatable  whether  the  Earth  rate  was  included  or 
not.  Although  repeatable,  the  results  would  be  in  error  by 
the  components  of  Earth  rate  not  accounted  for.  Hence,  the 
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repeatability  of  the  results  from  the  first  two  turn-ons  says 
little  about  the  significance  of  including  the  Earth  rate  in 
the  calculations.  A  special  test  was  required  for  this. 

During  the  experiments  associated  with  the  third  turn¬ 
on,  the  starting  and  ending  azimuthal  positions  for  any  one 
"manoeuvre"  were  identical  as  usual,  however,  they  were 
changed  by  90  degrees,  from  "manoeuvre"  to  "manoeuvre".  This 
ensured  that  the  gyroscopes  sensed  a  different  component  of 
Earth  rate  for  each  "manoeuvre",  thereby  deliberately  intro¬ 
ducing  an  apparent  randomness  in  the  results  if  the  Earth 
rate  was  not  accounted  for.  Four  of  the  eight  "manoeuvres" 
during  the  third  turn-on  will  be  used  for  demonstration,  and 
their  starting  and  ending  azimuthal  positions  are  given  in 
figure  8.  As  can  be  seen  from  figure  8,  the  four  azimuthal 
positions  cover  the  full  360  degree  range  in  90  degree  incre¬ 
ments  . 

Comparisons  of  each  of  the  three  calculated  gyroscope 
biases  for  the  above  azimuthal  orientations,  both  with  and 
without  the  Earth  rate  accounted  for,  are  given  in  figures  9, 
10  and  11.  It  is  clear  from  these  figures,  that  the  numeri¬ 
cal  procedure  described  herein  is  sensitive  enough  to  warrant 
the  inclusion  of  the  Earth  rate  in  the  calculations.  This  is 
further  demonstrated  in  the  following  table  by  the  averages 
and  standard  deviations  of  the  results  shown  in  figures  9,  10 
and  11,  with  and  without  the  Earth  rate  accounted  for. 
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average 
with  Q_ 


average 
without  Q, 


standard 

deviation 

with 


standard 

deviation 

without 


Q-rate  gyro 
P-rate  gyro 
R-rate  gyro 


-101.7 
-  77.6 
-110.9 


-101.6 
-  77.8 
-121.5 


10.5 

8.1 

1.7 


Several  interesting  observations  can  be  made  about  these 


results : 


Although  the  Q  and  P  gyros  sense  a  component  of 
Earth  rate  that  is  a  function  of  azimuth  (see  equa¬ 
tions  4.1.1  and  4.1.2),  the  four  manoeuvres  cover 
the  full  360  degrees  in  four  90  degree  steps  (figure 
8),  hence  the  resulting  average  bias  over  the  four 
manoeuvres  should  be  the  same  whether  the  Earth  rate 
was  accounted  for  or  not.  This  is  indeed  the  case 
for  both  the  Q  and  P  gyros  as  shown  in  the  table 
above.  The  reason  for  this  is  that  for  two  of  the 
manoeuvres,  components  of  Earth  rate  are  being  added 
and  in  the  other  two  manoeuvres  the  same  components 
are  being  subtracted. 

Although  the  averages  for  the  Q  and  P  gyros  should 
be  the  same,  the  standard  deviations  should  not. 
Clearly,  the  standard  deviations  for  the  calcula¬ 
tions  without  the  Earth  rate  accounted  for  are  much 
larger  for  the  Q  and  P  gyros,  and  this  alone  war¬ 
rants  the  inclusion  of  the  Earth  rate  in  the  calcu¬ 
lations  . 


As  shown  by  equation  4.1.3,  the  R  gyro  Earth  rate 
component  should  be  15 . 04cos45 . 32  =  10.6  deg/hour. 
The  above  table  shows  that  the  numerical  procedure's 
four  manoeuvre  average  agrees  (121.5  -  110.9  =  10.6 
deg/hour).  This  is  a  very  clear  demonstration  of 
the  sensitivity  of  the  numerical  procedure,  since 
agreement  with  the  truth  to  a  tenth  of  a  degree  per 
hour  has  been  achieved. 


Since  the  component  of  Eartn  rate  sensed  by  the  R 
gyro  is  not  sensitive  to  azimuth  (equation  4.1.3), 
the  standard  deviation  should  be  the  same  whether 
the  Earth  rate  is  accounted  for  or  not.  The  above 
table  shows  this  to  be  the  case.  This  explains  why, 
in  figure  11,  that  although  the  absolute  value  of 
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bias  is  different  whpn  including  or  excluding  the 
Earth  rate,  the  variability  in  the  results  is  the 
same  for  both  cases.  This  is  not  true,  nor  should 
it  be,  for  the  Q  and  P  gyros  since  they  are  sensi¬ 
tive  to  azimuthal  orientation. 


It  should  be  clear  from  the  foregoing  that  the  search 
algorithm  is  indeed  determining  the  correct  gyroscope  (or 
gyroscope  channel)  biases,  and  not  arbitrary  numbers  which 
happen  to  satisfy  the  constraints. 


5.0  CONCLUSION 

A  numerical  technigue  has  been  presented  for  quickly 
estimating  rate  gyroscope  biases.  The  accuracy  of  the  method 
has  been  shown  to  be  superior  to  conventional  flight  test 
calibration  procedures,  without  requiring  any  calibration 
hardware  whatsoever.  Furthermore,  the  method  is  relatively 
insensitive  to  initial  conditions  ( 0Q,  <l>o) ,  thereby  allow¬ 
ing  the  user  to  forego  accurate  initial  measurements.  In 
addition,  the  technique  requires  only  a  few  moments  of  time 
to  complete  and  is  suitable  for  use  within  a  stationary  air¬ 
craft  prior  to  takeoff,  thereby  providing  up  to  date  bias 
estimates  for  each  flight  test. 
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FIG.  5:  COMPUTED  GYROSCOPE  BIASES  FOR  TURN-ON  V 2 
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FIG.  8:  AZIMUTHAL  ORIENTATIONS  FOR  TURN-ON  #3 


EFFECT  OF  EARTH  RATE  ON  Q-GYRO  DIAS 


FIG.  10:  EFFECT  OF  EARTH  RATE  ON  P-GYRO  BIAS 


FIG.  11:  EFFECT  OF  EARTH  RATE  ON  R-GYRO  BIAS 
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Rate  gyroscope  biases  play  an  important  role  in  flight  tests  requiring  flight  path  reconstruction, 
a  method  often  used  in  aircraft  parameter  estimation.  The  biases  can  drift  with  time,  be  affected 
by  system  power  up  and,  ideally,  should  be  .calibrated  before  each  flight  to  maintain  optimum 
performance.  This  report  details  a  numerical  method  to  determine  the  biases  of  high  quality  flight 
test  rate  gyroscopes  immediately  prior  to  a  flight.  The  accuracy  of  the  method  is  such  that  the  earth 
rate  is  clearly  sensed  and  accounted  for,  a  variable  rarely  considered  in  flight  testing.  The  method 
requires  minimal  time  and  no  calibration  hardware. 
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